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Abstract 

We study physical situations when one or two "guest" arbitrarily- 
charged particles are immersed in the bulk of a classical electrolyte 
modelled by a Coulomb gas of ± unit point-like charges, the whole 
system being in thermal equilibrium. The models are treated as two- 
dimensional with logarithmic pairwise interactions among charged 
constituents; the (dimensionless) inverse temperature (5 is considered 
to be smaller than 2 in order to ensure the stability of the electrolyte 
against the collapse of positive-negative pairs of charges. Based on 
recent progress in the integrable (l+l)-dimensional sine-Gordon the- 
ory, exact formulas are derived for the chemical potential of one guest 
charge and for the asymptotic large-distance behavior of the effec- 
tive interaction between two guest charges. The exact results imply, 
under certain circumstances, anomalous effects such as an effective at- 
traction (repulsion) between like-charged (oppositely-charged) guest 
particles and the charge inversion in the electrolyte vicinity of a highly- 
charged guest particle. The adequacy of the concept of renormalized 
charge is confirmed in the whole stability region of inverse tempera- 
tures and the related saturation phenomenon is revised. 
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1 Introduction 



This paper deals with physical situations when one or two "guest" charges, 
say arbitrarily charged colloidal particles with a hard core of radius a, are 
immersed in a classical electrolyte modelled by an infinite Coulomb gas of 
positive/negative unit charges. In order to obtain explicit results we con- 
sider the point-like limit of the guest charges, i.e. a/A — ► where A is a 
characteristic correlation length of electrolyte species; the obtained results 
are not expected to be applicable to large-sized colloids. If the charges of the 
guest particles are sufficiently large, anomalous counterintuitive phenomena 
emerge in the system PQ. 

One of such phenomena is the appearance, under some circumstances, 
of an effective (i.e., mediated by the electrolyte) attraction between like- 
charged colloids. While the traditional DLVO theory [21 E] always predicts 
an effective repulsion between two like-charged colloids jU E] , experimental 
measurements (HI E] and numerical simulations (HI E] provide evidence for 
attraction, especially within confined geometries (close to a dielectric wall 
or between two glass plates) but also in the bulk of the electrolyte [10J. It 
was argued that the effective attraction of two like-charged colloids in the 
presence of a single wall can arise also from a non-equilibrium hydrodynamic 
effect PHH2]- 

Another interesting effect is the overcharging, or the charge inversion, of a 
highly charged colloid • This effect occurs when the number of electrolyte 
counterions in the vicinity of the colloidal surface becomes so high that the 
colloidal bare charge is locally overcompensated. The charge inversion has 
been observed experimentally by electrophoresis ^3] and in simulations jTH] • 
Its theoretical explanation is based on Wigner- crystal theories [To! IT7] . 

The latter effect is related to the concept of renormalized charge fHJ E3 
I2TH |2*H |22 I2H] ■ The true electric potential far from the colloid immersed in 
a weakly-coupled electrolyte is supposed to exhibit the Debye-Hiickel form, 
but with a renormalized-charge prefactor which is different from the bare 
charge of the colloid. An important feature, which occurs in the framework 
of the nonlinear Poisson-Boltzmann equation, is that the renormalized charge 
saturates monotonically at some finite value when the colloidal bare charge 
goes to infinity [221 12H]- Monte- Carlo simulations of a salt-free colloidal cell 
model [23] indicate the existence of a maximum in the plot of the renormal- 
ized charge versus the bare colloidal charge. Tellez and Trizac [23 considered 
the possibility of a more general phenomenon of potential saturation. 

A theoretical elucidation of anomalous phenomena requires to go beyond 
mean-field approximations by incorporating electrostatic correlations among 
electrolyte particles. Heuristic phenomenological approaches applied so far 
are based on plausible, but not rigorously justified, arguments. Some exactly 
solvable models are needed. The best candidates are two-dimensional (2D) 



2 



Coulomb systems with logarithmic pairwise interactions among the charged 
constituents. The 2D Coulomb gas of ± unit point-like charges is stable 
against the collapse of positive-negative pairs of charges at high enough tem- 
peratures, namely for /3 < 2 where j3 is the (dimensionless) inverse temper- 
ature or coupling constant. The collapse point (3 — 2, at which the collapse 
starts to occur, is equivalent to the free-fermion point of the Thirring rep- 
resentation of the 2D Coulomb gas \U\ ; although the free energy and 
the particle density diverge, the truncated Ursell correlation functions are 
finite at this point. In a recent work [21], we have solved exactly the 2D 
problem of one colloid immersed in the Coulomb gas just at the free-fermion 
point. An explicit form of the induced electric potential as a function of the 
bare colloidal charge was derived at every point of the space. Based on this 
exact result, the concept of renormalized charge was shown to fail in this 
strong-coupling regime. On the other hand, the anticipated phenomenon of 
the electric potential saturation was confirmed at the free-fermion point. 

Our present aim is to extend the exact treatment of the guest-charge (s) 
problem to the whole Coulomb-gas stability region of inverse temperatures 
< (3 < 2. A first important step towards this aim has already been done by 
solving exactly the equilibrium statistical mechanics of the 2D Coulomb gas 
in the stability regime (the bulk thermodynamics, special cases of the surface 
thermodynamics and the large-distance behavior of the two-body correlation 
functions) via an equivalence with the integrable 2D Euclidean sine-Gordon 
theory; for a short review, see ref. As is shown in this paper, the 

problem of one (two) guest charge (s) immersed in the Coulomb plasma is 
related to the evaluation of one-point (two-point) expectation values of the 
exponential field in the sine-Gordon theory. Based on recent progress in 
the latter topic, we derive explicit formulas for the chemical potential of 
one guest charge and for the asymptotic large-distance behavior of the effec- 
tive interaction between two guest charges. The exact results imply, under 
some circumstances, an effective attraction (repulsion) between like-charged 
(oppositely-charged) guest particles and the charge inversion in the elec- 
trolyte around a highly- charged guest particle. The adequacy of the concept 
of renormalized charge is confirmed in the whole stability region < (5 < 2. 
The related saturation phenomenon is revised. 

The paper is organized as follows. Basic facts about the bulk properties 
of the 2D symmetric Coulomb g clS £1X6 summarized in Section 2. Section 3 
deals with the problem of one guest charge in the electrolyte. The effective 
interaction between two guest charges immersed in the electrolyte is studied 
within a form-factor method for the equivalent sine-Gordon model in Section 
4. Based on the exact results of Section 4, the concept of renormalized 
charge and the related saturation phenomenon are tested in Section 5. A 
brief recapitulation and some concluding remarks are given in Section 6. 
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2 Basic facts about the 2D Coulomb gas 



We consider an infinite 2D plane A of points r G R 2 , filled with a homoge- 
neous medium of dielectric constant = 1. The electrostatic potential v at a 
point r, induced by a unit charge at the origin 0, is given by the 2D Poisson 
equation 

Av(r) = -2tt6(t) (2.1) 

The solution of this equation, subject to the boundary condition Vt>(r) — > 
as |r| — > oo, reads 

v(r) = -\n(^pj, reR 2 (2.2) 

The free length constant tq will be set for simplicity to unity. This definition 
of the Coulomb potential in 2D maintains many generic properties (e.g., 
sum rules |30j) of "real" 3D Coulomb fluids with the interaction potential 
v(t) = l/|r|, r G -R 3 . 

The symmetric Coulomb gas consists of two species of point-like parti- 
cles with opposite unit charges qj G {+1,-1}; to simplify the notation, the 
elementary charge e is set to unity. The bulk properties of the system in 
thermodynamic equilibrium are usually treated within the grand canonical 
ensemble. The ensemble is characterized by the (dimensionless) inverse tem- 
perature P, which plays the role of the coupling constant, and by the couple 
of equivalent particle fugacities z + = z_ = z. Since the length scale r in 
(|2.2j) was set to unity, the true dimension of z is [length]" 24 -^/ 2 ). The grand 
partition function of the plasma is defined by 



00 1 r N \ 
= E AT m j I a II [ d2r i z 9j\ ex P { 'I 3 E wMlrj - r k 

JV + ,7V_=0 7V + !7V - ! JA j=l I j<k 



(2.3) 

where (iV_) is the number of positively (negatively) charged particles 
and = + A^_. The system is stable against the collapse of positive- 
negative pairs of unit point-like charges provided that the corresponding 
Boltzmann factor r _/3 is integrable at short distances in 2D, i.e. for (3 < 2. 
In what follows, we shall restrict ourselves to this stability region of coupling 
constants. 

To introduce the averaged many-particle densities, we denote by (•••)$ 
the standard thermal average. At the one-particle level, one considers the 
number density of particles of one sign 



(EW( r - r i)) > ? = ±1 ( 2 - 4 ) 

Due to the charge symmetry and space homogeneity, n + = n_ = n/2 where 
n is the total density of particles. At the two-particle level, one considers 
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the two-body number densities 



n qq ,{r,r 



q,q' = ±l (2.5) 



which are translationally invariant, n, 



VI 



Mr, r 



n, 



QQ 



(|r-r'|). The two-body 
densities decouple at asymptotically large distance onto the product of the 
corresponding one-body densities, lim^oo n qq >{r) = n q n q >. It is therefore 
natural to introduce the Ursell functions, U qq i{r) = n qq i(r) — n q n q /, which go 
to as r — > oo. It is also useful to consider the pair distribution functions 
9 qq > (r) = n qq ,{r)/(n q n q ,). 

The short- distance behavior of the two-body densities is dominated by 
the Boltzmann factor of the corresponding pair Coulomb potential [HU E2] . 
In particular, the pair distribution functions behave like 



9qq'{r) 
C qq ' 



C nn >r M as r -> 
/3(/i- + ft - q ,) 



■qq 

exp 



(2.6) 
(2.7) 



Here, the excess (i.e. over ideal) chemical potential of the Coulomb-gas 
species is defined by 



cxp 



/4 X ) = -> 
' n q 



q = ±l 



(2.8) 



and fi'g with arbitrarily-valued Q represents an extended definition of the 
excess chemical potential: /ig is the reversible work which has to be done 
in order to bring a particle of charge Q (in units of the elementary charge 
e) from infinity into the bulk interior of the considered Coulomb gas. In the 
case of oppositely-charged particles, formula (|2.fij) reduces to 



9-\ 




as r 







(2.9) 



According to Eq. (J2TJ), -A/(2vr) is the inverse operator of the Coulomb 
potential v(r). The grand partition function of the 2D Coulomb gas (j2.3j) 
can be thus turned via the Hubbard-Stratonovich transformation (see e.g. 
ref. (231) i n t° 

^_JV<PeMS(z)} (21Q) 



with 



S(z) 



dV 



16n 



JV<f> exp[-S(0)] 



(j)A(j) + 2z cos(60) 



I 
4 



;2.ii) 



being the 2D Euclidean action of the sine-Gordon model. Here, <p(r) is a real 
scalar field and / T><p denotes the functional integration over this field. The 
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fugacity z is renormalized by the (diverging) self-energy term exp[/3u(0)/2], 
without changing the ^-notation. The one- and two-body densities of the 
charged particles in the plasma are expressible as averages over the sine- 
Gordon action as follows 

n q = z 9 (eW), rv(|r - r'|) = z q z q> (e^e^^) (2.12) 

With regard to Eq. (l2~5J) . it holds 

exp(-frif) = (e^>, q = ±l (2.13) 

The short- distance behavior ()2.9|) is equivalent to 

(e^We-^M) ^ | r _ r '|-^ 2 as | r _ r '| ^ o (2.14) 

Under this conformal normalization of the exponential field, the divergent 
self-energy factor (which renormalizes z) disappears from statistical rela- 
tions calculated within the sine-Gordon representation. The short- distance 
formula ()2.14j) is the special case of a more general relation 

(e^W^W) ^ ( e i(a+a '^)|r - r'| 4aa ' as |r - r'| -»• (2.15) 

valid in a restricted region of the parameters a and a' such that the one-point 
average (e l(a+ct '^) be finite. 

The sine-Gordon model (|2.11|) is integrable jM]. Its particle spectrum 
consists of one soliton-antisoliton pair (S, S) with equal masses M and of 
S — S bound states ("breathers") {Bf, j = 1,2, ... < l/£} whose number at 
a given b 2 depends on the inverse of the parameter 

The mass of the -breather is given by 

rrij = 2Msin (^3) ( 2 - 17 ) 

and this breather disappears from the spectrum just when rrij = 2M. The 
breathers exist only in a subinterval of the stability region < b 2 < 1/2 
(0 < (3 < 2) of the point-like Coulomb gas. The lightest i^-breather, usually 
called the elementary one, has the mass 

m 1 = 2Msin^y^J (2.18) 

and disappears from the particle spectrum just at the free-fermion point 
b 2 = 1/2 {(3 = 2). The soliton-antisoliton pair is present in the spectrum 
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up to the Kosterlitz-Thouless transition point b 2 = 1 (f3 = 4) at which the 
sine-Gordon theory ceases to be massive. 

The (dimensionless) specific grand potential u of the 2D Euclidean sine- 
Gordon model, defined by 



uj = lim — — -lnS 

|AHoo |A| 



(2.19) 



was found in ref. [HE] by using the Thermodynamic Bethe ansatz: 



— u 



8 sin«) 



(2.20) 



Under the conformal normalization of the exponential fields (J2.14)) . the re- 
lationship between the fugacity z and the soliton/antisoliton mass M was 
established in ref. 



7iT(l - b 2 ) 



M 



2T(£/2) 



2-2b 2 



(2.21) 



where T stands for the Gamma function. Note that the mass M has dimen- 
sion of an inverse length. As a consequence of Eqs. ()2.20|) and (j2.21|) . one 

ld(-uj) M 2 



(e ib ^) = - 



2 dz 8z(l - b 2 ) 



te 



Relations ()2.21|1 and ()2.22|1 . together with the equality 

<e^> 



n 
2~z 



(2.22) 



(2.23) 



determine explicitly the density-fugacity relationship and consequently the 
complete thermodynamics of the 2D Coulomb gas in the stability region j^Z] . 



3 One guest charge in the electrolyte 

Let us consider a point-like particle of charge Q with Q being an arbitrarily 
valued real number; when Q is interpreted as the valence it has to be an 
integer. The charge is put into the bulk interior of the 2D electrolyte, say 
at the origin 0. The electrostatic potential induced by the guest charge at a 
point r G R 2 is equal to — Qln|r|. Its effect on the constant species fugacities 
z± is the following: z q — > ^'(r) = ^Irl^ 9 . The excess chemical potential of 
the guest particle is given by 
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where S[z^(r)] represents an obvious functional generalization of the defi- 
nition (|2.3|) of the grand partition function with position-dependent particle 
fugacities. Performing the Hubbard-Stratonovich transformation, the sine- 
Gordon representation of 5[^W(r)] takes the standard form of Eq. (j2.10|) 
with the action 

SM(z) = - I d 2 r ( — + z\v\ AQb2 ^ + z\r\-^ b2 e~ ib A (3.2) 
J a KlQn J 

We first shift the scalar field 0(r) — > 0'(r) = 0(r) — i4Q£» ln|r | . Using subse- 
quently the Poisson equation Aln|r| = 27r<5(r) and integrating by parts with 
vanishing boundary contributions at |r| — > oo, one gets 

exp (-(3^) = (e iW > (3.3) 

where the average is taken with the usual sine-Gordon action S(z) given by 
Eq. (|2.11j) . When Q — ±1, one recovers the previous result (|2.13|) derived 
for the plasma constituents. Note the obvious symmetry /Zq = /i^g. 

An exact formula for the expectation value of the exponential field (e ia ^), 
where the sine-Gordon parameter b 2 lies inside the stability region < 
b 2 < 1/2 and a is a free real parameter, was conjectured by Lukyanov and 
Zamolodchikov in ref. [38J. In terms of our notation a = Qb [see Eq. ()3.3j) ]. 
their formula reads 

Q 2 b 2 

^ exp [I h {Q)\ , \Q\<^ (3.4) 



with 

w)=rf 

J t 

The formula was "guessed" on the base of three exactly solvable cases of the 
sine-Gordon theory: the semi-classical limit b 2 — > 0, the free-fermion point 
b 2 = 1/2 and the special value of a = b, see Eq. ()2.22|) . The validity 
of the formula was supported later by a "reflection" relationship with the 
imaginary Liouville theory j3U| > a numerical study of the sine-Gordon model 
in finite volume j^T] and a variational perturbation theory |22|. Other checks, 
provided by the Coulomb-gas representation, are presented in the next two 
paragraphs. 

The integral (Q is finite provided that \Q\ < l/(26 2 ); at \Q\ = l/(26 2 ), 
the integrated function behaves like 1/t for t — > oo what causes the loga- 
rithmic divergence of the integral. In the Coulomb-gas picture, the interac- 
tion Boltzmann factor of the Q-charge with an opposite unit plasma charge 
(counterion) at distance r, r~^\ is integrable at small r in 2D if and only if 



(e iW ) 



irzT(l - b 2 ) 



sinh 2 (2Q& 2 £) 



2 sinh(6 2 t) sinht cosh[(l - b 2 )t] 



-2Q 



2 b 2 e~ 2t 



(3.5) 
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(3\Q\ < 2. The stability region for /j,q therefore is expected to be \Q\ < 2/(3; 
there is a collapse at \Q\ =2/(3 characterized by /j,q — > — oo. With regard 
to the relationship (3 = 4b 2 we conclude that the pair of Eqs. ()3.4|) and ()3.5|) 
passes the guest-charge collapse test. 

The way in which (e 1 *^) diverges as \Q\ approaches the collapse value 
l/(26 2 ) is another check provided by the Coulomb-gas representation. To 
show this fact, let us attach a hard core of radius a around the guest Q- 
charge. The effect of the Q-charge on the species fugacities is the following: 



z g \r 



PQq, 



|r| — er), where 8 is the Heaviside step function. The 



procedure analogous to that outlined between Eqs. ()3.1|) and ()3.3|) leads to 



exp 



where the average is taken with the action 
S a (z) = S(z) + z [ d 2 r 



gi^(r) _|_ e -ib<t>(r) 



(3.6) 



(3.7) 



We expand (e 1< ^ f "^°' ) ) (J in the lowest <r-order around the sine-Gordon action 
S(z): 



(e iQbm ) a = (e lW > - z 



iQbcj>\ 



d 2 r 



r<(T 



iQ6</.(0) e i6</»(r)^ + ^ e iQ6<^»(O) e -i60(r) 



+ ■ 



(3i 



Due to the symmetry /Xq(ct) = ^^(a), it is sufficient to consider the case 
Q > 0. Applying in Eq. ()3.8)1 the short- distance formula (j2.15j) . the leading 
cr-correction is obtained in the form 



3 i<W(0)\ 



)a ~ {e iQH ) - nz(e 



a 



2-4Qb 2 



1 - 2Qb 2 



as a — > 



(3.9) 



Close to the collapse value of the guest charge, i.e. when Q = l/(2b 2 ) — e 
with e — ■> + , one expands 



a 2-AQb 2 a 4b 2 e 



1 - 2Qb 2 2b 2 e 2b 2 e 



+ 21na + 0(e) 



(3.10) 



Since the regularized (e 1 ^ 6 ^^)^ has to be finite at the collapse Q = l/(2b 2 ) 
it follows from Eqs. (J3.9)) and (J3.1U)) that the singular behavior 



lim ( exp 



1TZ 

( exp 



2b 2 t 



(3.11) 



must take place. It is shown in the Appendix that the singular behavior 
(I3.11J) is indeed reproduced by the conjectured Eqs. ()3.4|) and (|3.5p . 
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4 Two guest charges in the electrolyte 

Let us put two point-like particles into the bulk of the Coulomb plasma, the 
one with the charge Qi at the point ri and the other with the charge Q 2 at the 
point r 2 . The electrostatic potential induced by these two charges at a point 
r G R 2 is equal to — Qiln|r — ri| — Q 2 ln|r — r 2 |. The constant species fugacities 
z± are thus modified as follows: z q — ► zP^(r) = z\r — ri|^ l9 |r — r 2 |^ 29 . The 
excess chemical potential of the guest 1,2-charg GS clS db whole is given by 



cxp 



_ ^/ i QiQ 2 l ri ' r ^ -I r l- r 2| — ^73 — (4--U 



Under the Hubbard-Stratonovich transformation, the sine-Gordon represen- 
tation of S[z^ 2 )(r)] takes the standard functional form of Eq. (j2.10j) with the 
action 

S^(z) = - /d 2 rf— ^A^ + ^r-rxl^V-r^V^ 

J A \167T 

+z\v - ri r 4Qlb > - r2 |- 4 ^ 2 e-^ (4.2) 

Shifting the scalar field <p(r) — > <f)'(r) = <p(r) — i4Qi61n|r — ri| — i4Q 2 Mn|r — r 2 |, 
applying the Poisson equations Aln|r — rj| = 27r<5(r — rj) (j = 1,2) and 
integrating by parts, one arrives at 



exp 



(3Hq iQ M - r 2 )| = (e^^e"**^) (4.3) 



where the average is taken with the usual sine-Gordon action (j2.11|) . The 
effective interaction energy between the guest 1,2-charges is defined by 

iWh - r 2 |) = MQ lQ2 (|ri - r 2 |) - fa - ^ 2 (4.4) 
With respect to the relation (jH.Hjl . it holds 

/ e iQl^(ri) e iQ 2 fe0(r 2 )\ 

exp [-/3E QlQ2 (| ri - r 2 |)] = ^ QlH) ^ lQ2H) (4-5) 

At asymptotically large distance |ri — r 2 | — > 00, the two-point correlator 
^ e iQib0(ri) e iQ 2 b0(r 2 )^ decouples onto the product (e iQl ^) ( e iQ2 ^) and so the 
interaction energy goes to zero as it should be. From Eq. (|4.5j) one then gets 

/ e iQife0(ri) e iQ 2 fe0(r 2 )\ 

- PE QlQ2 (\ ri - r 2 |) ~ ^q^^qm) - 1 * ki - r 2 | ^ 00 (4.6) 

This means that the asymptotic large-distance behavior of the effective in- 
teraction energy between the guest particles is related to the large-distance 
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behavior of the corresponding two-point correlation function of exponential 
fields associated with the 2D sine-Gordon theory. 

For the 2D Euclidean sine-Gordon model, like for any integrable 2D the- 
ory, the two-point correlation function of local operators O a (a is a free 
parameter) can be formally expressed as an infinite convergent series over 
multi-particle intermediate states |43j . 



(Ojr)OAr')) = V - V I °° " F a {B u ■ ■ ■ ,9 




ejv—«l 



7N)ei-e N 



N 

r'| 5^m ej cosh0j ] (4.7) 



where e indexes the particles [say e = +(— ) for a soliton (antisoliton) and 
e = j for a £>j-breather] and the rapidity 9 G (— oo, oo) parametrizes the 
energy and the momentum of the corresponding particle. The form factors 

F a (9 1 ,...,9 N ) ei ... eN = (0|a(O)|Z £1 (^),...,Z £jv (^)) (4.8) 
eN -^F a/ (9 N ,...,9 1 ) = (Z eN (0 N ),...,Z ei (0 1 )\O a ,(0)\0) (4.9) 

are the matrix elements of the operator at the origin, between an iV-particle 
in-state (being a linear superposition of free one-particle states \Z e {9))) and 
the vacuum. The first N = term of the series expansion ()4.7j) corresponds 
to the decoupling {O a ){O a i). 

The form-factor representation (|4.7jl is particularly useful in the limit 
|r — r'| — ► oo where the dominant contribution to the truncated correla- 
tion function (O a (r)O a i(r')) — (O a ){O a i) comes from a multi-particle state 
with the minimum value of the total particle mass Y,f=i m ej , at the point 
of vanishing rapidities. As was already mentioned, the lightest particle in 
the stability region < b 2 < 1/2 is the elementary breather B\. For this 
particle, the one-particle form factors F a {9)\ and 1 F a i{9) = F a >(9)i of the 
exponential operator O a (r) = exp [ia0(r)] were calculated in refs. 



where 



F a {9)i = (Ole^ISi^)) = -i<e^)V^A Sm . (7 ^°{ 6) (4.10) 

A = -sin(7r0cosf^exp{- - — ) (4.11) 
7r y 2 J [ Jo f sint J 

and ^ is defined in Eq. ()2.16jl . Since the form factor (|4.1U|) does not depend 
on the rapidity, the integration over 9 in (|4.7j) can be done explicitly by using 
the relation 

00 d 9 / tt \ 1 / 2 

_ e -r mi coshe = Ko f mir \ ^ exp(-m 1 r) as r -> oo 

^oo 2 \2rriir J 

(4.12) 
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where K Q is the modified Bessel function of second kind [46] . On the base of 
the asymptotic equivalence ([4.6)1 . the large-distance behavior of the effective 
interaction energy is finally found in the form 

PE QlQ2 (r) ~ [Qi][Q 2 }\ (^) 1/2 exp(-m 1 r), r -> oo (4.13) 

Here, the symbol [Q] stands for the ratio 

_ sin(^g/(4-/?)) 
iWJ " sin(^/(4-^) l4 ' i4j 

Using the thermodynamic formulas derived at the end of Section 2, the mass 
rrii (which plays the role of the inverse charge-charge correlation length of 
the Coulomb-gas particles) is expressible as 



mi — K 



sin(7r/?/(4-/3)) 
n (3/ (4-/3) 



1/2 



(4.15) 



where k = y/2ir(3n denotes the inverse Debye length. The /3-dependence of 
the parameter A, defined by Eq. (|4.11j) . reads 

A = — sin cos — — — exp < — / : — > (4.16) 



2(4-/3) J { Jo vrsint 

An interesting feature of the result 1)4.13)1 is that the effective interaction 
energy factorizes into the product of separate charge contributions [Q] from 
each of the guest particles. 

We would like to emphasize that the asymptotic formula ()4.13)) was de- 
rived for the guest particles of point-like nature. Its rigorous validity is there- 
fore restricted to such guest charges which do not collapse with an opposite 
unit counterion from the electrolyte, i.e. to the values \Qi\, IQ2I < 2//3. Since 
the function [Q] is analytic at every real Q, the stability border \Q\ = 2/(3 
does not represent an exceptional point at which a singularity emerges [like 
it was in the case of the excess chemical potential determined by Eqs. (J3.3)) - 
(13. 5|) ]. The explanation of this important fact follows from the definition 
(14 .4|) of the effective interaction energy: if one of the guest charges passes or 
is beyond its collapse value, say \Q\\ > 2/j3, both excess chemical potentials 
yUg 1 Q 2 (|ri — r 2 |) and /^gj tend to —00 in such a way that their difference 
is expected to keep a finite value. We therefore suggest that the formula 
()4.13|) remains valid for arbitrary real values of Qi and Q2, and corresponds 
to the limit of a small hard-core radius a (such that niio << 1) around 
the guest charges. In what follows, we shall refer to this conjecture as "the 
regularization hypothesis" . 



12 



The interaction energy £?Q 1 Q 2 (r) is repulsive (attractive) at asymptoti- 
cally large distance r when it goes to zero from above (below). When the 
amplitudes of the guest charges are the same, i.e. \Qi\ = [Qzl, the inter- 
action energy (j4.13|) exhibits the vacuum-type behavior: it is repulsive for 
Qi = Q2 and, since [-Q] = —[Q], attractive for Q\ = — Q2- The situa- 
tion is more complex when the amplitudes of the guest charges differ from 
one another. Let us analyze the plot of [Q] as the function of (say posi- 
tive) bare charge Q for a fixed value of (3 from the stability range (0, 2). In 
the interval Q G (0, (2/(3) — (1/2)), [Q] increases monotonically from to 
its maximum at the end-point of this interval. In the subsequent interval 
Q E ((2/(3) - (1/2), (4/(3) - 1), which contains the stability border Q = 2/(3, 
[Q] is the decreasing function of Q but still keeps the positive sign of Q. 
Such behavior means physically that by increasing the bare charge of one of 
the guest particles the effective interaction energy weakens which is a coun- 
terintuitive phenomenon. The function [Q] changes the sign when Q passes 
the value (4/(3) — 1; under the assumption of validity of the regularization 
hypothesis, this indicates an effective change of the sign of the bare charge 
Q. Appropriate combinations of the Q±, (^-charges in the factorized relation 
()4.13|) can therefore lead to an effective attraction (repulsion) between like- 
charged (oppositely-charged) guest particles. The described scenario repeats 
itself when increasing Q due to the periodicity relation [Q] = [Q + (8/(3) — 2\. 

2.5 

2 
1.5 

1 

O 0.5 


-0.5 
-1 
-1.5 

1 2 3 4 5 6 

Q 

Fig. 1. The function [Q] versus the bare charge Q for various values of (3. 
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The above discussed plot of [Q] versus the (positive) bare charge Q is 
presented graphically for various values of the electrolyte coupling constant 
(3 in Fig. 1. The solid-line parts of the plots correspond to the interval 
< Q < 2/(3 for which the formula (j4.13j) is valid rigorously, the dashed- 
line parts of the plots correspond to Q > 2/(3 when the formula (j4.13|) is 
applicable under the assumption of the regularization hypothesis. 



5 Renormalized charge 

Let us put one point-like particle of charge Q at the origin and look for 
the evoked density profiles n q (r) of the electrolyte species q — ±1. One can 
formally follow the procedure for two guest particles (one of which has the 
charge q = ±1 of the electrolyte species) outlined at the beginning of Section 
4, to obtain 

/ e iQb^(0) e iqb<j,{r)\ 

n ^ r )= z J^q^ > ^ = ±1 (5.1) 

Note that in the special case Q = the spatially-homogeneous relation 
()2.23j) is reproduced. At asymptotically large distance r from the guest Q- 
charge, the straightforward application of the form-factor method explained 
in Section 4 provides the result 

n q( r ) ~ n q j 1 - [o\[Q]^ (2^7) exp(-mir)j, r -> 00 (5.2) 

where the symbols [q], [Q] are defined by Eq. (|4.14|) and the parameters mi, 
A by the respective relations ()4.15|) and ()4.16|) . 

The charge density p of electrolyte unit charges is defined as p(r) = 
n + (r) — n_(r). Since [q] — q for q — ±1, one finds that 

/ 7T \ l ' 2 

p(r) ~ —n[Q]\ exp(— mir), r — ► 00 (5.3) 

\2m 1 r/ 

The average electrostatic potential if) induced by the guest Q-charge is 
related to the charge- density profile through the 2D Poisson equation, 

A^(r) = -27rp(r) (5.4) 

Inserting here the asymptotic formula (|5.3|) and considering the circularly 
symmetric Laplacian A = r _1 d r (rd r ), the dimensionless electric potential 
(3ip can be shown to behave at large distance r from the Q-charge as follows 

mTJ ^ A V2mT/ exp(-m 1 r), r -> 00 (5.5) 
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In the Debye-Hiickel limit (3 — > and for finite Q, it holds [Q] ~ Q, mi ~ k 
and A ~ (3. Eq. (|5.5|) thus reduces to the well-known result (see e.g. refs. 

P4>Du(r) ~ 7; — exp(-Acr), r ^ oo (5.6) 



,2m", 

The asymptotic behaviors ()5.5|) and ()5.6|) . considered respectively in terms of 
the dimensionless combinations m\r and kt, exhibit the same type of the fall- 
off. This fact confirms the adequacy of the concept of renormalized charge 
|18| [23] in the stability weak-coupling regime of the Coulomb gas. Eq. ()5.5j) 
is consistent with Eq. ()5.6j) provided that one introduces the renormalized 
charge Q ren as follows 

PQren= ( ~ V [Q]A (5.7) 



This formula can be simplified to 



Q ren (P,Q) = A(P)sm(^^j (5.8) 
where the positive amplitude A((3) is given by 

-(4 - (3) sin —rr—^ ex P { I : — f ( 5 - 9 ) 



A(f3) 2 y rj \2(A-/3)J r [Jo Trsint 

As was mentioned above, because of the point-like nature of the guest 
Q-charge the rigorous validity of the result ()5.8|) is restricted to \Q\ < 2/(3. 
Changing Q from towards positive real values, the renormalized charge 
()5.8p increases monotonically up to its maximum A(j3) at the point Q = 
(2/(3) - (1/2). Increasing then Q from (2/(3) - (1/2) up to the stability bor- 
der 2/(3, Qren paradoxically decreases while still keeping the positive sign of 
the bare charge Q. This scenario resembles the one observed in the Groot's 
Monte-Carlo simulations of the salt-free (only counterions are present) col- 
loidal cell model |2U • 

Under the assumption of validity of the regularization hypothesis, one 
can further increase the value of Q beyond 2/(3 in the relation (|5.8|) . Qren 
changes its positive sign to the negative one at Q = (4//5) — 1. This change 
of the sign is closely related to the effect of charge inversion PUJ-JZj: si nce 
the total screening cloud of electrolyte particles must compensate exactly 
the bare charge Q of the guest particle (the electroneutrality sum rule), the 
fact that the charge density (|5.3|) goes to at asymptotically large distance r 
from above is the evidence of the charge inversion starting at some distance 
from the guest Q-charge. 

The renormalized charge Qren is a periodic function of Q. This is why 
going with Q — ► oo does not imply the saturation of Q re n at some finite 
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value. Instead, Q ren oscillates between the two finite ±A((3) extremes. With 
regard to Eq. (|5.5|) . the induced electrostatic potential exhibits the same 
type of the oscillatory behavior as Q — > oo what contradicts the idea of the 
monotonic electric potential saturation [25J. 

6 Conclusion 

Let us summarize briefly the crucial results of the present work. 

Section 3 deals with the case of one point-like particle of charge Q (with 
\Q\ < 2/(3) immersed in the bulk of the stable 2D Coulomb gas. Passing 
to the sine-Gordon representation, we were able to relate in Eq. ([3.3)1 the 
excess chemical potential /Xq of the guest charge to the expectation value of 
the exponential field. The explicit form of the latter quantity was conjec- 
tured by Lukyanov and Zamolodchikov [j3Ej, see Eqs. ()3.4j) and ([3.5)1 . and 
subsequently verified by various methods. Our Coulomb-gas formulation pro- 
vides two other checks of this conjecture: the guest-charge collapse test at 
| Q | = 2/(3 and the predicted singular behavior (|3.11|) of the exponential- field 
expectation close to the collapse value of Q = {2/(3) — e (e — > + ). 

The problem of two guest point-like particles, charged by Q\ and Q2 
and being at distance r from one another, is studied in Section 4. Using 
the form-factor method for the two-point correlation functions in the sine- 
Gordon formulation of the problem, we have derived the explicit formula 
()4.13|) for the effective interaction energy of the two guest charges Eq^qJ^t) 
at asymptotically large distance r — > 00. This formula is valid rigorously 
for \Qi\, IQ2I < 2/(3] in a subspace of this region of charge values we have 
noticed an anomalous weakening of the effective interaction when one of the 
guest charges is increasing. Since in the definition ([4.4)1 of Eq^^t) the 
divergences of chemical potentials (caused by the collapse of guest charge 
with electrolyte counterions) cancel with each other, we have suggested an 
extended validity of the formula ([4.13)1 for arbitrary values of Q± and Q2 (the 
regularization hypothesis). Then, under certain circumstances (especially, 
the inequality \Qi\ 7^ \Q 2 \ must hold and one of the guest charges has to 
be large enough), the result ()4.13[) implies an effective attraction (repulsion) 
between like-charged (oppositely-charged) guest particles. 

The adequacy of the concept of renormalized charge was confirmed in 
Section 5. This is not a surprise: in the whole stability interval of inverse 
temperatures < (3 < 2, the large-distance behavior of two-point correlators 
is determined by the form-factor of the same particle from the sine-Gordon 
spectrum, namely the i?i-breather with the lightest mass m\ playing the role 
of the inverse correlation length of electrolyte species. The large-distance 
behavior of the induced electric potential ([5.5)1 . considered in terms of the 
dimensionless combination mir, is therefore basically the same in the Debye- 
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Hiickel limit — > as well as at every point which belongs to the stability 
interval, up to the renormalized-charge prefactor. The renormalized charge 
Q r en, considered as a function of the (positive) bare charge Q, exhibits a 
maximum at Q = (2/0) — (1/2) which is in the range \Q\ < 2/0 of the 
rigorous validity of the formula ([5.8)1 . Under the assumption of validity of 
the regularization hypothesis, increasing Q can produce the effect of charge 
inversion. Going with Q — > oo in Eq. (|5.8j) does not imply the saturation of 
Qren at some finite value. Instead, Q ren oscillates between two finite extremes. 
The same property holds for the electric potential in the electrolyte region 
which contradicts the idea of the monotonic potential saturation [25J. 

The previous results obtained at the free-fermion (collapse) point (3 = 2 
|2*H] are different from the present ones concerning the stability interval 
< < 2 in the following aspects. Firstly, the concept of renormalized 
charge fails at — 2. Secondly, at — 2, when the bare charge Q — ► oo the 
induced electrostatic potential saturates monotonically at a finite value in 
each point of the electrolyte region. The reason for the fundamental differ- 
ences is obvious. The lightest .Bx-breather disappears from the particle spec- 
trum of the sine-Gordon model just at — 2, and the asymptotic behavior of 
two-point correlation functions at this coupling constant is governed by the 
soliton-antisoliton pair. Since mi — > 2M as — > 2, the particle mass in the 
exponential decay is a continuous function of at = 2. On the other side, 
the inverse-power-law asymptotic behavior, which is determined by the form- 
factor of the dominant particle(s) in the sine-Gordon spectrum, undertakes 
an abrupt modification when passing through the = 2 point. The basic 
qualitative features of the results obtained at the free-fermion point = 2 
are expected to be present also for such > 2 where the soliton-antisoliton 
pair exists. It is known [33] that the soliton-antisoliton pair disappears from 
the sine-Gordon particle spectrum (and the sine-Gordon theory ceases to be 
massive) at the point b 2 = 1 (0 = 4) which corresponds to the Kosterlitz- 
Thouless transition of infinite order from the conducting (fluid) phase to the 
insulating phase. We conclude that the 2D results obtained in the weak- 
coupling regime of the Coulomb gas < < 2 differ substantially from 
those in the strong-coupling regime 2 < < 4. 



Appendix 

According to the conjectured Eqs. ([3.4)1 and ([3.5)1 . it holds 



lim ( exp 

6^0+ ^ 



nzT{l - b 2 



46 2 (l-b 2 ) 



exp (h + h + h 



(a.i; 
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where 



h 
h 
h 



i dt 



sinh(t) 



1 



-2t 



o t [2sinh(6 2 t)cosh[(l - b 2 )t] 2b 2 
dt { 



sinh(t) 1 __ 2t 

T \2sinh(6 2 t)cosh[(l - b 2 )t] ~ ~ 2^° 



roc Af 

/ — exp(-46 2 et) = -C - ln(46 2 e) + 0(e) 
Ji t 



(A.2) 
(A.3) 
(A.4) 



and C is the Euler's constant. Using the same Eqs. ()3.4j) and ()3.5j) for 
expressing the expectation value of the exponential field on the rhs of Eq. 
(j3.11j) . one gets 



lim 

e^0+ 



where 



^exp 


i( 


Jb 


-eb)cj> 


) tiz r(i-6 2 ) 


( ex P 


i 


(± 

\2b 


-b) j 


\ " 2b 2 e T(b 2 ) 



exp (I[ + I' 2 -C- ln2) (A.5) 



1 Jo t \ sinh(t) 

00 dt j sinh[(l - 2b 2 )t] 



2^-2* 



l t 



sinh(t) 



-2(l-6 2 )e- 



■2t 



With the aid of the integral representations 



C 
ln2 



i dt 

o 7 
f°° dt / 

Jo TV 



1 - c" 



00 dt 
— < 
i t 



e '-e" 2i 



(A.6) 
(A.7) 

(A.8) 
(A.9) 



the argument of the exponential on the rhs of Eq. (jA.5|) can be written as 

dt f sinh[(l - 26 2 )t] 



I[ + I' 2 -C-\d2= I 
Jo 



t 



sinh(t) 



+ (2b 2 - l)e~ 2i 



(A.10) 



Finally, considering in Eq. (jA.5|) the integral representation of the logarithm 
of the Gamma function [JH] 

-(x-l)t _ 1 " 

(X - 1) + 



lnlYx) 



dt 



o t 



1-e- 



Re(x) > (A.ll) 



the proof of the desired formula (j3.11|) becomes accomplished. 



Acknowledgments 

I thank Bernard Jancovici for careful reading of the manuscript and useful 
comments. The support by Grant VEGA 2/3107/2003 is acknowledged. 



18 



References 

[1] Y. Levin, Rep. Prog. Phys. 65:1577 (2002). 

[2] B. V. Derjaguin and L. Landau, Acta Physicochim. URSS 14:633 (1941). 

[3] E. J. W. Verwey and J. Th. G. Overbeek, Theory of the Stability of 
Lyophobic Colloids (Elsevier, New York, 1948). 

[4] J. C. Crocker and D. G. Grier, Phys. Rev. Lett. 73:352 (1994). 

[5] J. Chakrabarti and H. L6wen, Phys. Rev. E 58:3400 (1998). 

[6] G. M. Kepler and S. Fraden, Phys. Rev. Lett. 73:356 (1994). 

[7] J. C. Crocker and D. G. Grier, Phys. Rev. Lett. 77:1897 (1996). 

[8] N. Gr0nbech- Jensen, K. M. Beardmore, and P. Pincus, Physica A 
261:74 (1998). 

[9] E. Allahyarov, I. D'Amico, and H. Lowen, Phys. Rev. E 60:3199 (1999). 
[10] N. Ma, S. M. Girvin, and R. Rajaraman, Phys. Rev. E 63:021402 (2001). 
[11] T. M. Squires and M. P. Brenner, Phys. Rev. Lett. 85:4976 (2000). 
[12] T. M. Squires, J. Fluid Mech. 443:403 (2001). 
[13] B. I. Shklovskii, Phys. Rev. E 60:5802 (1999). 

[14] H. W. Walker and S. B. Grant, Colloids and Surfaces A 119:229 (1996). 

[15] R. Messina, C. Holm, and K. Kremer, Phys. Rev. E 64:021405 (2001). 

[16] T. T. Nguyen, A. Y. Grosberg, and B. I. Shklovskii, J. Chem. Phys. 
113:1110 (2000). 

[17] R. Messina, C. Holm, and K. Kremer, Phys. Rev. Lett. 85:872 (2000). 
[18] G. S. Manning, J. Chem. Phys. 51:924 (1969). 

[19] S. Alexander, P. M. Chaikin, P. Grant, G. J. Morales, and P. Pincus, J. 
Chem. Phys. 80:5776 (1984). 

[20] H. L6wen, J. Chem. Phys. 100:6738 (1994). 

[21] A. Diehl, M. C. Barbosa, and Y. Levin, Europhys. Lett. 53:86 (2001). 
[22] E. Trizac, L. Bocquet, and M. Auboy, Phys. Rev. Lett. 89:248301 (2002). 



19 



[23] M. Auboy, E. Trizac, and L. Bocquet, J. Phys. A: Math. Gen. 36:5835 
(2003). 

[24] R. D. Groot, J. Chem. Phys. 95:9191 (1991). 

[25] G. Tellez and E. Trizac, Phys. Rev. E 68:061401 (2003). 

[26] F. Cornu and B. Jancovici, J. Stat. Phys. 49:33 (1987). 

[27] F. Cornu and B. Jancovici, J. Chem. Phys. 90:2444 (1989). 

[28] L. Samaj, Saturation of Electrostatic Potential: Exactly Solvable 2D 
Coulomb Models, cond-mat/0501442, to appear in J. Stat. Phys. 

[29] L. Samaj, J. Phys. A: Math. Gen. 36:5913 (2003). 

[30] Ph. A. Martin, Rev. Mod. Phys. 60:1075 (1988). 

[31] B. Jancovici, J. Stat. Phys. 17:357 (1977). 

[32] J. P. Hansen and P. Viot, J. Stat. Phys. 38:823 (1985). 

[33] P. Minnhagen, Rev. Mod. Phys. 59:1001 (1987). 

[34] A. Zamolodchikov and Al. Zamolodchikov, Ann. Phys. (N.Y.) 120:253 
(1979). 

[35] C. Destri and H. de Vega, Nucl. Phys. B 358:251 (1991). 

[36] Al. Zamolodchikov, Int. J. Mod. Phys. A 10:1125 (1995). 

[37] L. Samaj and I. Travenec, J. Stat. Phys. 101:713 (2000). 

[38] S. Lukyanov and A. Zamolodchikov, Nucl. Phys. 5 493:571 (1997). 

[39] D. Bernard and A. LeClair, Nucl. Phys. B 426:534 (1994); erratum 
Nucl. Phys. 5 498:619 (1997). 

[40] V. Fateev, S. Lukyanov, A. Zamolodchikov and Al. Zamolodchikov, 
Nucl. Phys. 5 516:652 (1998). 

[41] Z. Bajnok, L. Palla, G. Takacs and F. Wagner, Nucl. Phys. B 587:585 
(2000). 

[42] W. F. Lu, Phys. Lett. B 602:261 (2004). 

[43] F. A. Smirnov, Form-Factors in Completely Integrable Models of Quan- 
tum Field Theory (World Scientific, Singapore, 1992). 



20 



[44] S. Lukyanov, Mod. Phys. Lett. A 12:2543 (1997). 
[45] S. Lukyanov, Phys. Lett. 5 408:192 (1997). 

[46] I. S. Gradshteyn and I. M. Ryzhik, Table of Integrals, Series and Prod- 
ucts, 5th ed. (Academic Press, London, 1994). 



21 



